An Analytic Equation of State for Ising-like Models 
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Using an Environmentally Friendly Renormalization we derive, from an underlying field theory 
representation, a formal expression for the equation of state, y = f(x), that exhibits all desired 
asymptotic and analyticity properties in the three limits x —* 0, x — > oo and x — » — 1. The only 
necessary inputs are the Wilson functions 7a , 7 V and 7^,2 , associated with a renormalization of the 
transverse vertex functions. These Wilson functions exhibit a crossover between the Wilson-Fisher 
fixed point and the fixed point that controls the coexistence curve. Restricting to the case N = f , 
we derive a one-loop equation of state for 2 < d < 4 naturally parameterized by a ratio of non-linear 
scaling fields. For d — 3 we show that a non-parameterized analytic form can be deduced. Various 
asymptotic amplitudes are calculated directly from the equation of state in all three asymptotic 
limits of interest and comparison made with known results. By positing a scaling form for the 
equation of state inspired by the one-loop result, but adjusted to fit the known values of the critical 
exponents, we obtain better agreement with known asymptotic amplitudes. 

PACS numbers: 64.60.Ak, 



I. INTRODUCTION 

The universal equation of state for the Landau- 
Ginzburg- Wilson O(N) model remains a subject of great 
interest (see, for instance^ andii for recent reviews). Its 
calculation from first principles is much more difficult 
than the calculation of other universal quantities, such 
as critical exponents or amplitude ratios. The univer- 
sal equation of state for the scaling function exhibits 
crossover behavior between three distinct, asymptotic 
regimes - the critical region when approached from the 
critical isotherm, or when approached along the critical 
isochor, and the coexistence curve. Maintaining the cor- 
rect analyticity properties of the equation of state in all 
these three distinct regimes has not been possible within 
the confines of first principle calculations, such as from 
a field-theoretic, microscopic Hamiltonian. Rather, ap- 
peal has been made to a parameterized phcnomcnological 
scaling ansatai 3 - that exhibits the right asymptotics, the 
underlying microscopic theory then being used to fix var- 
ious free parameters that exist in the ansatz. 

In this paper we use a renormalization group methodol- 
ogy - Environmentally Friendly Renormalization- - based 
only on an underlying Landau- Ginzburg- Wilson Hamil- 
tonian and without the need for any phenomenological 
ansatz, to derive an equation of state that exhibits all 
desired analyticity properties in the three distinct asymp- 
totic regimes. 



II. THE EQUATION OF STATE 

As first noted by Widomi, in the critical region, the 
equation of state is a homogeneous function relating an 
external magnetic field H, the reduced temperature t, 



and the magnetization ip, which can be expressed by 

y = f{x) (i) 

where f(x) is universal, the scaling variables y and x are 
y = B s c H/tp s and x = B^t/p 1 /! 3 , and B c and B are 
non-universal amplitudes associated with the behavior 
on the critical isotherm t = and the coexistence curve 
t < 0, H = 0. 

The scaling function f(x) is normalized such that 
/(0) = I on the critical isotherm, and /(— 1) = on 
the coexistence curve. Several properties of the universal 
equation of state are known rigorously. For instance, it is 
known that f(x) has a regular Taylor expansion around 
the limit x = given by 
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n=l 



0^n 



(2) 



while in the limit x — > 00, by Griffith's analyticity, one 
has an expansion of the form 



/(*) = 



OO 



n=0 



00 -2n/3 



(3) 



In the limit x — > 00 a natural variable is z = bi(p/t@, 

where 61 = (— C 4 f /(C + ) ) and the C^ n are the ampli- 
tudes of the 2n-point correlation functions for T > T c , 
C + being the amplitude of the susceptibility. In terms of 
z the equation of state takes the form 



H = -{-C + /C+) 1/2 t 0S F(z) 



(4) 



where the universal scaling function F(z) for small z has 
an expansion of the form 



F{z) 



1 3 
z + -z 6 + 



E 



Tlr. 



J2.U-X 
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(5) 
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where r 2 = = 1 by choice of normalization. As JSJ) is 
an expansion in tp, the constants r 2n are related to the 2n- 
point correlation functions at <p = and hence are very 
natural observables to calculate in lattice simulations. In 
the limit z — » oo, -F(z) has an expansion of the form 

oo 

F(z)=z s J2 F k°z- k/0 (6) 

fc=0 

The universal scaling functions f(x) and F(z) are related 
via 

z- s F(z) = F^f{x) (7) 

with z — zqx~P, where zq is the universal zero of the 
equation of state in terms of the variable z. Hence, the 
expansion coefficients of the two functions can be related 
to find 

foo _ 2n+l-S r 2n+2 

Jn ~~ F °°(2n+1)! [ ' 

zpoo 

fn = ^0 n/ " 0) 

Thus, we see it is sufficient to know the expansion co- 
efficients of f(x) in the limits x — > and oo in order 
to calculate the asymptotic properties of F{z) and the 
interesting coefficients r 2n - 

Unlike the limits x — > and oo, near the coexistence 
curve, x — » — 1, there are no rigorous mathematical argu- 
ments as to the analyticity properties of f(x), although 
there do exist conjectures. For instance, for N > 1 ini^, 
based on an e-expansion analysis, it was conjectured that 
(l + x) has a double expansion in powers of y and y( d ~ 2 )/ 2 
of the form 

1 + x = ay + c 2 y 1 - e/2 + d iy 2 + d 2 y 2 ~ e/2 + ... (10) 

In three dimensions it predicts an expansion of (1 + x) in 
powers of y 1 ' 2 . 

Of course, the behavior in the vicinity of the coexis- 
tence curve depends on the value of N. For N = 1, 
the longitudinal correlation length remains finite away 
from the critical point on the coexistence curve, while 
for N > 1, the existence of Goldstone bosons leads to in- 
frared singularities. For N = 1, one may formally posit 
that in the vicinity of the coexistence curve 

oo 
n=l 

the integer powers being a reflection of the finite longitu- 
dinal correlation length. For Ising-like systems essential 
singularities are to be expected. Of course, these cannot 
be captured within the confines of an ansatz like 1)110. 
For N > 1 studies of the non-linear a model lead one to 
expect a leading behavior of the form 

f(x)~ Cf (l + x) 2 W-V (12) 



though, as mentioned, the nature of the corrections to 
this behavior is not well understood, although (|l(Jfl is 
one conjecture. In the l/N expansion there is some 
evidence 11 for logarithmic corrections of the form ln(l + 
x) in three dimensions. 

Early field theoretic calculations using the renormal- 
ization group with an £-expansioni4 foundered on the 
fact that they did not exhibit Griffiths analyticity in the 
large x limit. Irrespective of the expansion method - £-, 
fixed-dimension, l/N etc. - there will remain a funda- 
mental difficulty - that an expansion around a particular 
fixed point will not readily access the universal equation 
of state in the entire phase diagram, due to the presence 
of other fixed points that must be accessed. Simply put, 
the equation of state exhibits crossovers, and the nature 
of these crossovers depends on N. For N = 1 the theory 
is controlled by a "Gaussian" or mean-field fixed point 17 , 
wherein fluctuations are suppressed on the coexistence 
curve away from the critical point by the non-vanishing 
longitudinal mass. In distinction, for N > 1, the theory 
is dominated by the massless Goldstone excitations on 
the coexistence curve and the non-linear a model gives a 
good description 2 . 

The problem of incorporating Griffiths analyticity was 
solved by using a parameterized formulation in terms of 
new variables R and 9, related to t and (p via 

ip = m^R 13 m{9) 
t = R(l-6 2 ) 
H = hoR 0S h{9) 

where the two functions m(9) and h{9) are undetermined. 
In this new parameterization the scaling variable x and 
the scaling function f{x) are given by 

Most current field theoretic formulations for determin- 
ing the equation of state (se o 1 ' 11 for comprehensive re- 
views) rely on such formulations. The drawback is that 
the underlying microscopic theory is not used to deter- 
mine the functional form of m{9) and h{9), rather an 
ansatz is made as to the general functional form which 
depends on certain unknown parameters and then the 
underlying microscopic theory is used to fix these pa- 
rameters. The most common ansatz is that the func- 
tions are polynomials in 9. The coefficients of the pow- 
ers of 9 are then determined by calculating certain ob- 
servables independently from the underlying microscopic 
theory and then using the values of these observables to 
determine the coefficients. Various methodologies have 
been usediii*i£ using a variety of methods. For instance, 
an expansion of the effective potential for small ip for 
t > using e or fixed dimension expansions. Results 
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from Monte-Carlo simulations or high temperature ex- 
pansions have also been used. 

In this paper we apply a new method to describe the 
equation of state in a uniform manner. Using environ- 
mentally friendly renormalization 9 , we find a schema able 
to capture the crossover between the Wilson-Fisher fixed 
point and the fixed point that controls the coexistence 
curve. By integrating along curves of constant magneti- 
zation, we obtain the equation of state in the whole criti- 
cal region. Moreover, the equation of state that we obtain 
is parameterized in terms of the inverse of the transverse 
correlation length, a quantity well defined over the en- 
tire phase diagram. The representation we have found 
is valid for both large and small values of the scaling 
variables and satisfies Griffiths analyticity. Although our 
approach works for any N we here concentrate on the 
case N = 1. 



III. A RENORMALIZATION GROUP 
REPRESENTATION OF THE EQUATION OF 
STATE 

In this section we briefly outline the derivation of the 
equation of state for a theory described by the standard 
LGW Hamiltonian with O(N) symmetry 



must be imposed, where ft is an arbitrary renormalization 
scale. The renormalized parameters satisfy the differen- 
tial equations 



K[<p] 



with r = r c + tg, where r c is the value of r at the critical 
temperature T c and ts = A 2 ^ T Z Tc ' , A being the mi- 
croscopic scale. We denote a generic vertex function by 
t[ N \^\ , where the number of / and t subscripts indicates 
whether a longitudinal or a transverse propagator is to be 
attached to the vertex at the corresponding point. When 
all subscripts are either I or t we will use a single I or t, 
for example T 



(N,M) 



will be abbreviated T 



(N,M) 



Further- 



more, when there are no ip 2 insertions (i.e. M = 0) the 
second index will be left off e.g. indicates r^' ). 

Due to the Ward identities of the modeli&, it is suf- 



ficient to know only the T 



(N,M) 



as all the other vertex 



functions can be reconstructed from these. The equation 
of state is 



h = r 



(2) 



(16) 



A. Renormalization in terms of Non-Linear Scaling 

Fields 



Due to the existence of large fluctuations in the critical 
regime a renormalization of the microscopic bare param- 
eters of the form 



t(m, ft) 
A(ft) 
<p(n) 



Zj(K)t B (m) 
Z\{k)\ b 



(17) 
(18) 
(19) 



dtln) . . . . 

dA(re) 



dn 
dip(n) 



7a(k)A(k) 
1 



dn 



7 v (ft)^(ft) 



(20) 
(21) 
(22) 



where the Wilson functions associated with this coordi- 
nate transformation are 



7 V 2 («) 
7a(«) 
7 V («) 



— k— In Z 

dn 

K^-\nZ x 

dft 

ft— \nZ v 
dn 



■-p- 



(23) 
(24) 
(25) 



and the derivative is taken along an appropriately chosen 
curve in the phase diagram, which we here denote by 
c. Similarly, integration of the renormalization group 
equation for any multiplicatively renormalizable r^'^ 
yields 



-,(JV,M) 



(*(ft),A(ft),^(ft)) (26) 



To impose a specific, as opposed to abstract, coordinate 
transformation between bare and renormalized theory 
the renormalization constants Z v , Z v 2 and Z\ must be 
fixed. Here, we impose the explicitly magnetization de- 
pendent normalization conditions 



<Vr[ 2) (p, i(K, ft), A(k), ip(k),k) 

p 2 =0 

r< 2 '%t(K,K),A(K),¥'(K),K) 
,(4) 



1 



(27) 
(28) 



rp(0,i(ft,ft), \(k),<p(k),k) = A. (29) 

Note that in this case we impose the normalization con- 
ditions on the transverse correlation functions. These 
conditions serve to fix the three Z functions associated 
with ifB, ts and Xb while the condition 



ft 2 = t[ 2 ^ (0, i(ft, ft), A(ft), </?(k), ft), 



(30) 



serves as a gauge fixing condition that relates the slid- 
ing renormalization scale ft to the physical temperature 
t and the physical magnetization ip. Physically, k is a 
fiducial value of the non- linear scaling field m t , which is 
the inverse transverse correlation length. 

Besides m t , the other non-linear scaling field we use to 
parameterize our results is 



1 



(2) 



(31) 



lp 2 =0 
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which is an RG invariant. It represents the anisotropy in 
the masses of the longitudinal and transverse modes and 



is related to the stiffness constant p s — tp d p zT 



(2) 1 



via mz 



lp 2 =0 



y 3 . 

With a given renormalization prescription one may de- 
termine the equation of state in terms of the non-linear 
scaling fields m t and the transverse and longitudi- 

nal propagators that appear in all perturbative diagrams 
can be parameterized in terms of them. One of the main 
motivations for this reparametrization in terms of non- 
linear scaling fields is that it eliminates all tadpole dia- 
grams at higher loop order. However, what is required is 
the equation of state in terms of the linear scaling fields 
t and tp. One must therefore determine the coordinate 
transformation, t = t(mt,m v ), ip = tp(mt,m v ), between 
them. 



one. In this case, dt integrates to the temperature vari- 
able (T — T c (f)), where T c (ip) corresponds to that point 
on the coexistence curve where the magnetization is tp. 
As we wish to use as temperature variable t = (T — T c ), 
we may write {T-T c (tp)) = t+A, where A = (T c -T c (ip)) 
is the temperature shift that measures the distance be- 
tween the critical point and the point on the coexistence 
curve, T c (ip). In the integral, — Aj v 2dx/x, one may 
safely take the universal limit, kq — > oo. In this universal 
scaling limit, in terms of the linear scaling fields t and tp, 
the problem has only one scaling variable, x. In terms 
of the coordinates m t and m v , this manifests itself as 
a reduction to the single scaling variable, z — m t /m v . 
Passing to this variable, using lpTT)l. and integrating this 
along a curve of constant tp (taking into account the Ja- 
cobian factor) one finds 



B. Relating Non-linear and Linear Scaling Fields 

This may be done by specifying a particular curve, c, 
in the phase diagram along which we integrate the dif- 
ferential relation for the tranverse vertex functions t[ N ^ . 
For example, 



dr\ 2) = r t 



(32) 



for rj ' and T\ ' ' and the normalization conditions l|28|) 
and H30JI one may write 



dT 



(2) 



(2,1) 



r; 



(2-7*)« 



" ^K Q 7 ^ 



' Kdn 



(33) 



In the universal limit, where A — > oo—, the crossover to 
mean field theory is pushed off to infinity and the theory 
is then controlled in the limit n — > oo by the Wilson- 
Fisher fixed point. Hence, ji —> , where 7^ F is 
the Wilson function 7, at the Wilson-Fisher fixed point. 
Defining A 7i = (7, - 7l WF ) yields 



-7^2 



dt={2- 7v )e~ J »o A7 - 2 ^ l^- J ndn (34) 

To integrate (|3*1J) we need to fix some boundary condi- 
tion. The coexistence curve, where m t = 0, is a natural 



i4i(l + aj) = T[z) 



(35) 



where the scaling variable x = B 1 '^t/(p 1 '^, B being the 
non-universal amplitude introduced previously, and the 
universal scaling function T(z) is 



T(z) = f 

JQ 



2 - 7a + l v 



where 



can be integrated along a curve of constant (p to yield D(x) = exp 

(2) (2 1) 

dt = dT t /T t ' , where the right hand side is naturally 
written in the coordinate system (m t ,m v ). To integrate 
the renormalization group equation for the vertex func- 
tions it is most natural to use m t = k as the flow variable 
and hold m v constant. However, if we then wish to inte- 
grate (|32fl along a curve of constant ip we must include a 
Jacobian factor, 2/ (2 — 7^ + j v ), that takes into account 
the variation of m v along a constant <p curve. The re- 
lation between m v and tp is specified by (|31|l using the 

renormalization group equations for d p 2 T ^ and with 
the normalization conditions H27|) and 129|) . Using H2ti|) 



A7. 



A 7 a , A 7v \ ^ 



2 - 7A + 1v> 



(36) 



The quantity Ai is related to the universal zero of the 
equation of state. In terms of the amplitude B 



B 2 = 



3 K (4- d -^ /i i 



(37) 



Equation (|35|l determines the coordinate transforma- 
tion, t = t(mt,mtp) <p = Lp(jn t ,m v ), in the scaling limit 
where there is only one relevant scaling variable associ- 
ated with the linear scaling fields, x ~ t/ip 1 /^ , and one 
relevant scaling variable associated with the non-linear 
scaling fields, z = m t /m v . Hence, we determine the co- 
ordinate transformation x — x(z). Note that the geom- 
etry of this transformation has some unusual properties 
relative to the (T, tp) plane. The coexistence curve is 
mapped to the single point x = — 1, the critical isotherm 
is mapped to the single point x — 0, and the critical 
isochor, for t > 0, maps to the single point at infinity. 
However, the critical point itself - as the intersection of 
the coexistence curve, critical isotherm and critical iso- 
chor - maps to all these points. 

To proceed to the equation of state we use (|16fl and 

(2) 

the renormalization group equation for T t . Once again, 
passing to the variable z and including in the Jacobian 
factor, one obtains 



7 = (^) G[Z) 



(38) 



5 



where the universal scaling function Q{z) is given by 



Q(^Z^ = Z^e 2 "7A+T^ y g 2 — 7A+7V W 



(39) 



Introducing the scaling variable y = A 3 1 (H/ip d ), where 



A 3 is related to the non-universal amplitude B c via 

(3-1) 



Bt = A, 



A 
3^ 



(40) 



one sees that the universal equation of state is now of the 
form Q with 



f{x) = ±-Q{F- 1 {A l {l + x)) 

^3 



We may now ask what more can be said about A\ and 
A 3 , or whether they are simply non-universal parame- 
ters, being related to the non-universal amplitudes B c 
and B, that cannot be determined within the present 
formalism? A\, in fact, is related to the universal zero 
of the equation of state, and may be calculated in the 
following manner: We choose some arbitrary value of z, 
zq, and write J-{z) — > J-(zq) + F(x)dx; then choose 
z and zq in the asymptotic regime z, zo ^ oo z > zo, 
wherein, from Griffith's analyticity, we can write T{z) — 
z 1 ^ J2^Lo 3~™z~ 2n - In the limit z — > oo, depending on 
the value of j3, certain terms in the asymptotic expansion 
diverge. For instance, for the three dimensional Ising 
model, at one loop, (3 — 3/10, hence, only the n = 
and n — 1 terms diverge, while the contributions from 
n > 2 — > 0. Denoting the divergent part of the expan- 
sion as J- 00 ^); as zo is constant we may identify the non- 
constant term, Aix, in this limit with J-°°{z). Hence, we 
may identify A\ = lim Zo ^ 00 (JT(zo) — JT°°(zo)). In terms 
of the integrand, I{x), of the scaling function J 7 , we may 
write 



Note that this equation of state has been determined 
from a first principles calculation based on an underlying 
microscopic model. It is parameterized, but parameter- 
ized in a way that is completely determined by the un- 
derlying model. This is in distinction to standard para- 
metric represent ation o 1 1 1 1 , where, after imposing certain 
analyticity requirements, there is a large arbitrariness in 
determining the scaling functions h(6) and m(9). In fact, 
these functions may depend on an arbitrary number of 
parameters and, for each parameter, a universal quan- 
tity must be independently calculated in order to fix it. 
Also, the parameter z here has a much more transpar- 
ent and direct physical meaning than 9, being simply the 
ratio of the two fundamental non-linear scaling fields in 

(41) 

the problem - m t and m v - the transverse correlation 
length and the stiffness constant. As these quantities 
are well defined throughout the phase diagram this for- 
mulation has an added advantage relative to parametric, 
fixed-dimension expansions where the relevant non- linear 
scaling field used is the mass for T > T c and there are 
difficulties reaching the ordered phasei. 

In order to determine the expansion coefficients f° 
and as introduced in section 1, one requires the 
Taylor expansion of f(x) around x = and x = oo. 
In terms of our parametric representation, d n f{x)/dx n 
can be expressed using d/dx — (dz/dx)d/dz, where 
dz/dx — A\j{dT{z)jdz) hence, 



At 



(I(x) - I°°{x))dx 



(42) 



where I°°(x) is defined via T°° — J Q Z I°°{x)dx. A\ is 
clearly universal. To determine A3, we set the condition 
y = 1 on the critical isotherm t = 0. This corresponds 
to a particular value, z c , of z. Hence, A3 — Q(z c ). To 
determine z c note that from H36[) on the critical isotherm 
A\ = T(z c ). The inversion of this function allows for the 
identification of z c . Hence, we deduce A3 to be 



A 3 = G(T- 1 (A 1 )) 



(43) 



which is, once again, a universal function. 

So, given that A3 is determined from A±, and A\ is 
determined from T by subtracting off its divergent com- 
ponent as z — ► 00, we see that the simple ingredients 
that enter into a complete specification of the universal 
equation of state are the three Wilson functions - 7a, 7^2 
and 7^. These are the only quantities that need to be 
determined perturbatively (or otherwise). 



d n f(x) 
dx n 




(44) 



which need to be evaluated at the points of interest z — z c 
(x = 0), z = 00 (x = 00) and z = (x = — 1). For 
instance, taking the limit x — > 00 in (|36|l and (|39|l . and 
using the fact that the Wilson functions 7, approach their 
values at the Wilson-Fisher fixed point so Aji — > 0, one 
finds 



T(z) -> 7Z 



V/3 



Hence, 



/(*) - ^7™ 7 ^ 7 



+ 0(x^) 



(45) 
(46) 



(47) 



from which, using (|43|) we may identify the expansion 
coefficient 



./o 



(Ai/7) 7 



(48) 



which is related to the universal amplitude ratio R x via 
f§° = R^ 1 . Using the expression for C + in our formu- 



lation, C+ = K- (2 ~^ 



7 



and equations l|3*7jl and jlity. 



one can verify that l|48|) is identical to the expression 
R x = (C+fB^iB/Bc) 5 - 1 . Similarly using the expres- 
sion for C4" in our formulation 



(49) 
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one may determine the universal amplitude ratio R% 
-C+B 2 /{C+f to be 



Ri 



7 



7— di/ 



(50) 



With these two amplitude ratios in hand two-scale factor 
universality implies that any other may be determined. 



IV. ONE-LOOP RESULTS 

The advantages of the present formulation can be best 
illustrated by considering a concrete example. We will 
consider the universal equation of state in the one loop 
approximation, as at this level it is still possible to obtain 
analytic or "quasi" -analytic results. 

We begin with the values of the Wilson functions to 
one loop. The running dimensionless coupling A satisfies 



dX(z) 
dz 



eX + c d X (z) (1 + ^)— + 



1,4=2 , (N-l) 



9 



(51) 

where c d = 3(4 - d)T((4 - d)/2)/2{An) d / 2 . Taking the 
initial condition X(zq) = A, in the limit zq — > 00, X — » 00 
one arrives at the universal separatrix solution^ 



X(z)=[c d [(l + ^) * 



On the separatrix 

7a = (4-d) 



(52) 

(53) 

(54) 
(55) 



In the limit z — > 00, the Wilson-Fisher fixed point is ap- 
proached and ji — > 7^ F with, at one loop, 7a = (4 — d) 
and 7^2 = (4 - d)(N + 2)/{N + 8). On the contrary, 
in the limit z — > the strong-coupling fixed point is ap- 
proached and 7, — > 7^ c . For N > 1 the Goldstone bosons 
dominate and 7a = 7^2 = (4 — d). For N = 1 however, 
this fixed point is mean-field like as fluctuations are sup- 
pressed and 7i — > 0. 

In the limit z — ► 00, the Wilson functions can be ex- 
panded as power series in z~ 2 for any iV 




7i(«) = 7f F + > ai(n)z" 



(56) 



Hence, the universal scaling functions T and Q can also 
be written as power series in z~ 2 . This is true in a dia- 
grammatic expansion to all orders not just at one loop. 
The limit z — > is more complicated. In this limit, 



7, — > 7? c but the nature of the corrections is not ob- 
vious. At the one- loop level, from l|53|) . one can see that 
the leading corrections to the strong coupling fixed point 
values will be z( 4 ~ d )/ 2 . 

For N 7^ 1 analytic progress is difficult. However, for 
N = 1 these expressions simplify greatly yielding 



7A = (4-d) (l + l 



(57) 
(58) 



With these expressions one can explicitly calculate the 
scaling functions Q{z) and T{z) 



(59) 




(60) 



To determine the constant A\ we take the z — > 00 limit 
of (|60|l . identify the divergent part with A±x and the con- 
stant remainder with Aim^ Note that the case d = 2 is 
problematic. This is an artefact of the one-loop approx- 
imation where, in particular, in d = 2 one has /3 = 0, 
which is clearly unphysical. The case d — 4 is also inac- 
cessible due to the fact that we took the universal limit 
A — > 00. It may be recovered by returning to l|51|l and 
integrating it and not taking this limit. In this case, 
as expected, logarithms appear. So, equations (|59|l and 
(|60|l are valid for 2 < d < 4. A3 can be determined as a 
function of the undetermined A\ using l|43|l . Explicitly, 
taking the z — * 00 limit Ai can be identified to be 



d + 2 V d-2 



2(4-d) 
3(d-2) 



(61) 



from which one determines z c = 2 -1 / 2 which, interest- 
ingly, is dimension independent, though, once again, this 
is a one-loop artefact. With z c in hand one determines 
A 3 to be 



A3 = ^ 



(4-d) 

If d + 2 \ < d - 2 ) 



2 \2(d-2) 



(62) 



Thus, the one-loop equation of state for N = 1 in d- 
dimensions is 



2< d - 2 ) n n 

—z 2 (2 + (d-2)z 2 )— (63) 



(d+ 2)< d - 2 > 
(2z 2 - 1)(1 + (d- 2)z 2 /2)^^ 



(64) 



Griffiths analyticity can be seen quite simply from these 
expressions. For x — > 00 one has z — > 00 and the right 
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hand side of takes the form z 1 /' 3 X)^Lo( c ™/ z2 ™)' 

where /? = 3(d - 2)/2(d + 2). Substituting into 
one then has the expansion J2J), where 7 = 6/(d+ 2). 
Similarly, by expanding (|64f> around z = z c one obtains 
the expansion (0). Finally, in the vicinity of the coexis- 
tence curve, from B64JI we see that (1 + x) can be written 
as a power series expansion in z 2 and, hence, we verify 
pi[l. Once again, within the approximation used one 
would not expect essential singularities to appear. To 
summarize: (|63[1 and 1|64|) have been determined from 
first principles from the underlying microscopic theory 
and obey all analyticity and other properties required in 
the different asymptotic regimes. 

^From these expressions one may calculate analytically 
the coefficients /° (d) and /£° (d) as functions of d, as well 
as derived quantities, such as the ri ni and universal am- 
plitude ratios, like R x and Rj . In the same way, we may 
examine the behavior in the vicinity of the coexistence 
curve. First, we determine z as a power series in (1 + x), 
and then we replace the result in the expansion of y for 
z small. 

In Table [I] we show the coefficients fi(d) — /5(d) and 
in Table El/x(«0 - f§( d )- The coeffients / °°(d) - 
for arbitrary <i appear in Table UTTI Our results for some 
important quantities like r2 n {d) appear in Table ITVl Ob- 
viously with analytic expressions in hand it is straight- 
forward to generate other coefficients. 



suit by substituting d — 4 — e in (|63|l and l|64|) , obtaining 

2 

V = 7-^^-^(2(1 + z 2 )-ez 2 )(^> (65) 
(6 - e)^l 

x = (2z 2 -l)(l + (2-e)z 2 /2) T?T ^ 7T (66) 
Expanding in powers of e to O(e) one finds 



y = 2z A 1 + - In 



£ , 2(1 + z 2 ) 



(67) 



1+ar = 2z 2 + -(2z 2 -l)ln(l + z 2 ) (68) 
o 

Inverting (|68|l in powers of e and substituting into l|67|) 
one finds 

y = l + .T + -(l + x)ln^-^--xln^-^ (69) 

which is the well known result in terms of normalized 
variables. 



A. The Case d = 3 

For d = 3 things are even more transparent, as the 
parameter z can be totally eliminated, thereby ending 
with a direct, unparameterized relation between x and 
H/ip s . For d = 3 the scaling functions T[z) and C?(z) are 



TABLE I: Values of /° for the d- dimensional Ising universality 
class. 



3X2 3^3(d-2) (2 + d)( 



4+2 d 8d jg 

(rf-4) (2 + d)^^ 2 3 ( d ~V 



(1/3) x 2 



-3+7 



(rf-4) (2 + rf)^ 



7u-ia a 

(1/27) x 2 3d- 6 (d- 



7) (d- 



4) {2 + d)i^A (8 + d) 



(d-16) (ri-4) (5 + d) (2 rf-11) 
20 (d-4) 5 (2 + d) 

405X 2 3(d-2) ( 2 +d)3(d-2) 



Hz) = ^ + jj(, 2 + 2)§(2z 2 -l) (70) 



= (z 4 + 2z 2 ) 

while Ai and A3 are given by 

2§3 
M = — 



A, 



(71) 



(72) 



With these values for the amplitudes A\ and A 3 the 
equation of state in terms of the variable z is 



(73) 



TABLE II: Values of for the d-dimensional Ising univer- 
sality class. 



1± 

IS 



3-23=2 (2 + d)^ 



3 . 2^3=2- -3 ■ 2^" + 5 • 23=2 



d— 4 

(2+d) 



5 (2 + d) — 



d— 1 , , 4— d / 3d+4 2 + 4d 

64^ 7^=^(2 + d)3=5 81 • 2^=r + 47 . 2^2- 



(2+d)3 



d 2 (34-d) 2 + 4 d 

45 ■ 323=2 _ 35 . 2~ d + 9 • 2=i=r d 



One may also recover the well known e-expansion re- 



which can be simply solved (the positive square root is 
required) and substituted into Q7Uf) to find 



23 a: = 1 




(74) 



which can be seen to satisfy y — 1 at x — and y = 
at x = —1. For large x, y — * (2 8 / 5 /5).t 6 / 5 , i.e. y ~ x 1 
as required by Griffith's analyticity. It is valid for both 
t > and t < 0. 

In the second column of Table we see the numeri- 
cal values of the coefficients /° and f^°, as well as some 
important derived quantities, such as the r2„. All these 
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TABLE III: Values of /£° for the d-dimensional Ising univer- 
sality class. 



ft 

1' X 

J2 

/oc 
3 

roc 
Ji 



41^2 



d —4 

-2 



14-5 d 8-2d \ jqr^ 

2 e -- id (d-2) 



3 (d-4) d~T~ 

1 - - ' \ 2+d ( 2+d\*=2 

1^2 , 



2 6-3d (rf _ 2)3d-6 



3x2 ^ ^ 2 w ((i _ 2) ^ - |±d 



K(d-3) 



2~3^2" 



(d-4)(4d-13) 



I |z^ J .( 2-^3=3- (d- 2) 1 ,^_] 



2 +<* / 2+d\ 3=2 



(17 -5d) (lid -32) 



2(8d-35j 50-15d d^4 
2 d-2 (d-2) <<-a (|±|)3=5 
BB X 



45 2 6 -3 d (d-2) 3d- 



2d\ 2 + d 



/ 6d + 8 

176553 x d + 



27(2419d- 11973) d 2 2~ 



189 (1520 + 13d 4 ) 2~^~ 




FIG. 1: The scaling function f(x) of the three-dimensional 
Ising model. The dashed line was taken from reference-. 



V. FITTED EXPONENTS 



TABLE IV: Values of , Fq° and some r 2 „ coefficients for 
the d-dimensional Ising universality class. 



R4 


3 (d-2) 

e / 5, 4 8-2 d\ 2+d 
_l_( 2 3 + ^f (d~2)^e) + 


F§° 


d-1 4-d 
42-d (d-2) 2 (§±4) d -2 

6 (d-3) 

/ 14-5 d 8-2 d\ ,4_ rf 
9 2 6-3d ( d _2)3d-6 


re 


|(4-d) (d-2) 


ra 


f (d-4) (4d-13) 


rio 


f (4-d) (5d-17) (lid -32) 


ri2 


385 x (665 x 2 4+ ^A _ 6539 x 2 6+ ^ d+ 

2+4 d 2+4 d d . \ 

11973 x 2^ d 2 - 2419 x 2^r d 3 + 91 x 323=5 rf 4 j 



values are as in good an agreement with known values as 
one might expect from a one-loop calculation, in those 
cases where a comparison can be made, with one appar- 
ent exception: the value of r 8 is about 2 — 3 times bigger 
than the majority of estimates, which are in the range 
2.18 — 2.7. However, a Monte Carlo simulation of Kim 
and Landau^ led to a value of r§ which was much larger 
than other estimates. However, this estimate goes down 
when we fit with the exact exponents in the next section. 
Our value of rio = —17.5 - is in the expected region of 
previous calculations, which give estimates in the region 

—4 25, though the errors associated with many of 

these estimates are large. Our estimate of r\2 for d = 3 
would seem to be new. Obviously an important advan- 
tage of the present approach is the facility with which the 
T2n can be calculated for even very high n. In the figure 
we see a comparison between our one-loop equation of 
state and that obtained by the HT method. 



Although one of the chief advantages of the present 
methodology is the fact that it is completely self- 
contained, in that there are no parameters to be fixed 
by appealing to exogenous information, as in standard 
parametric approaches, it is possible adapt the present 
method to utilize information that is available, such as 
precise estimates for the critical exponents and amplitude 
ratios. To illustrate this here, we once again consider the 
case N = 1. We take the one-loop "crossover" function 



(1 



") to give the exact form of the crossover and fit 



its asymptotic value to the best estimates for the critical 
exponents. Hence, we take 



7A 



(4-d) 



7<p2 



(2-1/iz) 



(1 



(75) 



The crossover form for rj is of course a pure supposition 
given that the form (1 + z~ 2 ) is derived from a one-loop 
calculation. The form of (|75|l is such as to guarantee a 
crossover to the known asymptotic behavior. The same 
procedure could be carried out using a two-loop calcula- 
tion. In this case the crossover functions for each Wilson 
function would be different, due to the fact that differ- 
ent diagrams contribute to them. Once again constants 
would be introduced to ensure a crossover in the limit 
z — > 00 to the correct exponent values for v and 77. With 
the ansatz (|75[) one finds 



Hz) 



is\^ 2/3(7-1) 







1 - 2/3 



7(1-2/3) 2 
21/(7 " 1) 



2,i ' 



g{z) = ^ ( 1 + _ 



(76) 
(77) 
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Examining the large z limit one determines the universal 
amplitude A\ to be 



Consequently, we determine = 2^(7 
and 



_ 2/3(7-1) / 2^(7-2/3) 
3 (7-2/3) ^27/3(1-2/3) 



(78) 

1)/ 7 (1 - 2/3) 



(79) 



Using this mechanism, in the third column of Table 
fVl we show the results for the three-dimensional Ising 
model. The values of the critical exponents that we have 
used are the best values reported in the literature 
7 = 1.2372, /3 = 0.3265 and v = 0.6301. 
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TABLE V: Numerical values of expansion coefficients for the 
three-dimensional Ising universality class: HT results taken 
from reference—, one loop(7?Gl) and fitted exponents (RGA) 
results. The resulting values for the coefficients rm were ob- 
tained from the relation J5J and using the values of /^°. 





HT 


RG1 


RGA 


TOO 

Jo 


0.6024(15) 


0.606 


0.596 


roc 
Jl 




0.696 


0.793 


roc 
72 




0.6 


0.613 


roc 

h 




0.223 


0.151 


/»oo 
74 




-0.066 


-0.084 


roo 

JS 




-0.0454 


0.00527 


re 


2.056(5) 


2.5 


1.938 


rs 


2.3(1) 


5.833 


2.505 


rio 


-13(4) 


-17.5 


-12.599 


ri2 




-192.5 


10.902 


/? 


1.0527(7) 


1.034 


1.05 


fS 


0.0446(4) 


0.029 


0.043 


fO 
J3 


-0.0254(7) 


-0.0034 


-0.0054 


f 

H 




0.0007 


0.0013 


Is 




-0.00019 


-0.0004 


ft 


0.9357(11) 


0.96 


0.939 


n 


0.08(7) 


0.048 


0.076 


fz 




-0.0112 


-0.023 


ft 




0.0049 


0.0138 


fi 




-0.0028 


-0.011 




7.81(2) 


6.892 


7.981 


F °° 


0.03382(15) 


0.0347 


0.0263 



being substantially different to the theoretical value of 
f{f = 0.6369— The only other notable change is that of 
/5 , which is related to ri2, where there is a sign change 
passing from the one-loop to the adjusted values. 



f (X) 



FIG. 2: The scaling function f(x) of the three-dimensional 
Ising model using fitted exponents. The dashed line was taken 
from reference-. Differences are not visible at this scale. 



TABLE VI: Numerical values of expansion coefficients for the 
two-dimensional Ising universality class; high precision (HT) 
results taken from reference— and results from using fitted 
exp onent s ( RGA ) 





HT 


RGA 


roc 

Jo 


0.14753 


0.1043 


rcc 
Jl 




0.1696 


rcc 

h 




0.2365 


J3 




0.2931 


r& 


3.67867(7) 


2.8571 


rs 


26.041(11) 


15.2381 


rio 


284.5(2.4) 


125.714 


ri2 


4.2(7) x 10 3 


1436.73 


fl 




1.1724 


fS 




0.1473 


fO 




-0.0164 


fx 




0.785 


fi 




0.295 


R4 


7.336774(10) 


9.7594 


^0 


5.92357(6) x 10" 5 


0.10067 



Most of the values found with fitted exponents sub- 
stantially improve the value compared with the one-loop 
approximation in those cases where a comparison can be 
made with HT expansions. This is particularly notable 
in r 6 , r 8 and rio . /g° and are a little puzzling in this 
respect. However, it is worth noting the experimental 
result reported in referenceii where /g° = 0.5917, this 



In the same way, we can also substitute the well- 
known exact values for the critical exponents in the two- 
dimensional case, i.e. 7 = 7/4, /3 = 1/8 and v — 1. Table 
121 shows the subsequent results. In this case, except for 
the value of /q°, the values found for r 2n show signifi- 
cant differences. This probably hints at the inadequacy 
of the ansatz for the crossover functions (|75|l . We have 
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not been able to find reported values for the coefficients 
/° and /£, some of which are reported in Table fVTl It 
would be interesting to be able to make the comparison. 

VI. CONCLUSIONS 

Using environmentally friendly renormalization we de- 
rived a formal expression for the equation of state for the 
O(N) model. This expression has the advantages that: 
i) it is derived from an underlying microscopic (field- 
theoretic) model; ii) requires as input for any calculation 
only the three crossover Wilson functions associated with 
a magnetization-dependent renormalization of the field, 
ip, the composite operator, ip 2 , and the coupling con- 
stant, A. In particular no experimental input is required 
on order to fix parameters; iii) it is parameterized by 
two non-linear scaling fields - the transverse mass, and 
an anisotropy parameter closely related to the stiffness 
constant; iv) it manifestly expresses all relevant, desired 
analyticity properties, both in the critical region and on 
the coexistence curve; v) universal coefficients associated 
with the expansion of the equation of state in any of the 
asymptotic regimes may be simply calculated thereby ob- 
taining coefficients that are presently unknown. 

After deriving one-loop expressions for general N, the 
formulation was then used to calculate a parameterized, 
analytic expression for the equation of state for N — 
1 for 2 < d < 4. For d = 3 it was shown how the 
parameterization could be dispensed with and a closed- 
form, non-parameterized expression for the equation of 
state derived. Various universal coefficients associated 



with the asymptotic regimes x — 0, x — oo and x = 
— 1 were derived for arbitrary d. For d = 3 these were 
compared with previous results. 

By taking the functional form for the Wilson functions 
at one-loop and introducing constants to fit to the best- 
known asymptotic values of the critical exponents, com- 
parison was made for d = 3 and d = 2 between our cal- 
culated expansion coefficients and those, where known, 
as derived using HT expansions, Monte Carlo, etc., in 
the different asymptotic regimes. In most cases the fit- 
ted values were substantially better than the one-loop 
values, agreement being better for d = 3 than d = 2. 

When continued to complex values of the external mag- 
netic field the universal equation of state should also cap- 
ture the Lee- Yang edge. There are additional points of 
non-analyticity in our expressions at z\ Y = — 2/(d — 2) 
or z\ Y = —v/f). These imaginary values of z are nat- 
urally associated with the Lee- Yang edge singularities. 
However, our equation of state is not yet optimized to 
include the associated crossover and we do not expect to 
obtain good estimates for the associated universal expo- 
nents or amplitudes. Our formulation can be adjusted 
to include this additional singularity, but as of yet we 
have not studied the effect of including this crossover. 
We hope to return to this in the future. 
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